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Abstract. We present the first results on the hard X-ray continuum image (up to 15 keV) of the supernova remnant Cas A 
measured with the EPIC cameras onboard XMM-Newton. The data indicate that the hard X-ray tail, observed previously, that 
■ extends to energies above 100 keV does not originate in localised regions, like the bright X-ray knots and filaments or the 

primary blast wave, but is spread over the whole remnant with a rather flat hardness ratio of the 8— 10 and 10—15 keV energy 
bands. This result does not support an interpretation of the hard X-radiation as synchrotron emission produced in the primary 
shock, in which case a limb brightened shell of hard X-ray emission close to the primary shock front is expected. In fact a weak 
rim of emission near the primary shock front is discernable in the hardest X-ray image but it contains only a few percent of the 
hard X-ray emissivity. The equivalent width of the Fe-K line blend varies by more than an order of magnitude over the remnant, 
it is hard to explain this as Fe-emission from the reverse shock heated ejecta given the ejecta temperature and the age of the 
remnant. The uniquely high wavelength-dispersive RGS- spectrometer has allowed, for the first time, to extract monochromatic 
images in several highly ionised element species with high spectral resolution. We present here a preliminary result on the 
measurement of the O VIII Ly-a and Ly-/3 brightness distribution and brightness ratios. The large observed decrease of the Ly- 
a/Ly-/3 ratio going from the N to the SE can be explained by small-scale (10 ") variations in the Nh column over the remnant 
and the potential presence of resonance scattering of the O VIII Ly-a photons in the limb brightened shell. 
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1 . Introduction tion spatially resolved spectroscopy of the remnant in the X-ray 

. . . domain is a powerful tool to study in detail the distribution and 

The young galactic supernova remnant (SNR) Cassiopeia A . . , . * , \ u i u . j i ^, 

,^ ax- .11 i i. i i , 1^/ ii physical properties of the (reverse) shock heated plasma. The 

(Cas A) is widely believed to be the result of the core collapse , . . ... , . . - . £ vata* at . 

v . J r bandwidth (up to 15 keV) of XMM-Newton and the presence 



of a massive star, probably an early type Wolf-Rayet star (Fesen £ .. , . .. . ' / /npC x ff 

i nr^k ^ a • i -n / -i of high wavelength-dispersive spectrometers (RGS) offers an 

et al. 1987^). Cas A is classified as an oxygen rich remnant . t . ( . . T ... 

1 — r-v, . , . (nC r v unique opportunity to address these questions. In this respect 



since optical spectroscopic observations (Chevalier & Kirshner 



XMM-Newton strongly complements Chandra, and this paper 



19791) show the supernova ejecta (in the form of fast moving . c ... ... u+ . u . v , 

\ r J v i gives a first highlight ol this capability 



knots) to contain mostly oxygen and oxygen burning products 
such as sulphur, argon and calcium. At all wavelengths Cas A 

has the appearance of a broken shell with a radius varying be- 2. Observations 

tween 1.6' to 2.5'. 



One of the outstanding problems in the study of Cas A 



A description of the instrument is given by Jansen et al. (2001 ). 



. ^ . . £ .i i- i i i v The data were obtained m July 2000 with a net exposure time of 

is the origin ol a recently discovered hard X-ray continuum _ . r r t 

. ^ . j. ^ . A „ Al A . , . A j. 30 ks. The telescope was pointed at the centre of the remnant 

m the spectrum of this remnant. From a theoretical point of , _ 09ho9moKB ^- oo ; o/ ^ . ^ n 



view both synchrotron emission from shock accelerated elec- 
trons and non-thermal Bremsstrahlung from electrons acceler- 
ated from the tail of the thermal emission could be possible ^ 1 

. . r» . . . t jj v i-i i The raw data were processed with the development ver- 

explanations for the hard continuum. In addition, high resolu- . _ , ^ ^ K ^ T « . A , . ^ ;« A „ N m 

sion of the XMM-Newton Science Analysis System(SAS). The 



(a = 23 h 23 m 25 s 5=58°48 / 20 // ) and the telescope roll angle 
was such that the RGS dispersion axis was aligned at 45° (NE) 
on the sky. 



Correspondence to: J.A.M.Bleeker@sron.nl RGS spectra were extracted by applying spatial filters to the 
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spectral image and the spectral order (m = —1) is selected ap- 
plying the appropriate pulseheight intervals to the CCD spec- 
tral camera. 

3. The hard X-ray data 

A high-energy tail has been observed in the X-ray spectrum of 
the supernova remnant Cassi opeia A by the Compton Gamma 
Ray Observatory (The et al. |1996|), BeppoSAX (Favataetal. 



1997) and the Rossi X-ray Timing Explorer (Allen et al. 1997) 



Previous hard X-ray im aging observations, for example using 
BeppoSAX (Vink et al. |1999|) indicated that the hard contin- 



uum radiation originated predominantly in the W region of the 
remnant, however this result was based on deconvolved images 
from a typically V resolution telescope with moderate effective 
area. The combination of large collecting area in the energy 
band 4.0 to 12.0 keV and angular resolution of a few arcsec- 
onds of XMM-Newton provides us with a unique opportunity 
to search for the distribution and origin of this hard "tail". 

3. 1. The X-ray spectrum above 4.0 keV 

Below 4.0 keV the observed spectrum is a complicated combi- 
nation of spectral lines, continuum emission and the effects of 
interstellar absorption but above 4.0 keV it is dominated by a 
smooth continuum and the Fe-K emission line complex. A joint 
spectral fit of the EPIC MOS and PN spectra was performed 
in the energy band 4.0—12.0 keV for the complete remnant. 
The model included a single Bremsstrahlung component plus a 
power law continuum with photon index fixed at the previously 
observed value (Allen et al. 1997) of 1.8 below a break energy 
of 16 keV. Gaussian line features were used to fit the obvious 
Fe and Ni line blends. The resulting fit is shown in Fig. |l[ The 
best fit Bremsstrahlung temperature is 2.56 ± 0.05 keV with 
emission measure 1.99 ± 0.04 x 10 59 cm -3 for a distance of 
3.4 kpc and the power law flux at 1 keV is 0.0416 ± 0.003 pho- 
tons cm -2 s _1 keV -1 . The Fe line energy is 6.603 ± 0.001 keV 
in the MOS and 6.623±0.002 keV in the PN. The errors quoted 
are statistical only and indicate the high quality of the data, the 
difference between the MOS and PN values may be attributed 
to systematic uncertainties which are at present of the order of 
10 eV. We estimated the effect of potential pile-up in the PN 
camera to be not larger than ~3 % for the brightest regions. 
Fig. [j] indeed shows that the MOS spectra, which do not suf- 
fer from pile-up, and the PN spectrum are fully compliant. The 
equivalent width of the Fe line is 1.02 keV for the MOS and 
1.11 keV for the PN. The Ni line energy is 7.740 ± 0.015 keV 
in the MOS and 7.783 ± 0.015 keV in the PN corresponding 
to Ni xxvi 7.687 keV and/or Ni xxvn 7.799 keV. The equiva- 
lent width of the Ni line was 0.130 keV for the MOS and 0.179 



keV for the PN. Allen et al. (1997) fitted the composite spec- 
trum of Cas A over a much broader energy range extending 
out to 100 keV using two Raymond-Smith collisional ionisa- 
tion equilibrium (CIE) thermal components with kTi=0.6 and 
kT2=2.9 keV to fit the soft band below 6 keV and a broken 
power law (i^=15.9 keV) with photon index 1.8 below and 
3.0 above E b with a flux of 0.038 photons cm -2 s _1 keV -1 at 




1 keV. Favata et al. (1997), from the BeppoSAX observations, 



channel energy (keV) 

Fig. 1. Joint spectral fit of the hard tail from the complete rem- 
nant, MOS 1 red, MOS 2 green and PN blue. The dashed and 
dotted line indicate the Bremsstrahlung and power law compo- 
nents contributing to the MOS continuum. 



obtained a continuum parameterized by two non-equilibrium 
ionisation (NEI) components with kTi=1.25 keV, kT 2 =3.8 keV 
and a power law with photon index 2.95 with a normalisation of 
0.74 photons cm -2 s _1 keV -1 at 1 keV. It is most important to 
note that, irrespective of the exact parameterization of these fits, 
they all share the fact that the X-ray flux above 8 keV becomes 
dominated by the power law component. More quantitatively, 
our spectral fit predicts a ratio between the power law and the 
Bremsstrahlung component in the MOS cameras of 0.29 in the 
interval 4-6 keV, 1.04 in the 8-10 keV band and 2.00 in the 
10—12 keV hardest band. We shall come back to this point in 
the next section on the X-ray images. 

3.2. X-ray images above 4.0 keV 

The MOS and PN images are practically identical. Fig. [| shows 
MOS images and Fig. || shows combined MOS and PN im- 
ages. The lefthand panel of Fig. || is the continuum map for 
4.07—6.07 keV an energy band in which no bright spectral lines 
are visible in either the EPIC MOS or PN spectra. The peak of 
the hard continuum occurs in the large knot to the W but there 
is also extended emission throughout the volume of the rem- 
nant. The outer rim (presumably just behind the blast wave) is 
visible but is very faint compared with the central volume. 

The righthand image of Fig. || is the Fe K line equivalent 
width derived from the narrow energy band 6.20—6.92 keV. 
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Fig. 2. Continuum 4.06-6.07 keV (left) and Fe K equivalent 
width (right) black keV to white 4 keV and above. 



The continuum used to estimate the equivalent width was es- 
timated by interpolation from the bands 4.06—6.07 keV and 
8.10—15.0 keV. The Fe line is very bright relative to the con- 
tinuum on the outer extremities of the SE knots, i.e. region 
marked P, but is rather weaker in the W where the continuum is 
bright, i.e. region marked L. The lefthand panel of Fig. ||is the 




Fig. 3. Continuum 8.10—15.0 keV (left) and hardness ratio 
(10.0-15.0 MOS+PN)/(8.10-10.0 MOS+PN) 0.15 black to 
0.45 white. 

hard continuum above the Fe K line, 8.10—12.0 keV, which is 
very similar to the softer continuum image in Fig. [| In order 
to try and locate where the hard X-ray tail originates we cal- 
culated the hardness ratio (10.0- 15.0)/(8. 10- 10.0) using both 
the MOS and PN cameras to give the best possible statistics. 
The result is shown in the righthand panel of Fig. ||. The con- 
tour superimposed in green is taken from the hard continuum 
map 8.10—15.0 keV and contains 37 % of the total counts in 
this energy band, the hardest region in the W contains less than 
3 % of the total. The remaining 63 % of the hard flux is spread 
out over the rest of the remnant. There are significant changes 
in hardness over the remnant but all the bright features within 
the green contour have a remarkably similar spectrum above 
8.0 keV. 

From the spectral fits described in the previous section, it 
is clear that the photon flux observed in the continuum image 
in Fig. ||is dominated by the hard tail (power law) component. 
XMM-Newton clearly detects the hard X-ray tail from the rem- 
nant but the hard X-ray image and the hardness ratio indicate 
that this flux does not predominate in a few localised regions, 
but pervades the whole remnant in a distribution similar to the 
softer thermal components. 



Table 1. Results of joint spectral fits using a single 
Bremsstrahlung and gaussian line 4.0— 12.0 keV. Listed are the 
temperature kT, line centroid Et, equivalent width EW and 
line sigma width a. The first values for Eg, EW and a are for 
the MOS cameras and the second from the PN. 



Parameter 


all 


L 


P 


kT (keV) 


3.53±0.02 


3.41±0.14 


3.15±0.34 


E t (keV) 


6.60 


6.60 


6.67 




6.62 


6.61 


6.68 


EW (keV) 


0.95±0.02 


0.57±0.08 


6.19±0.56 




1.01±0.02 


0.43±0.18 


5.28±0.53 


a (keV) 


0.083 


0.079 


0.047 




0.078 


0.042 


0.031 



3.3. The Fe K emission 

Simpler spectral fits using a single Bremsstrahlung and one 
gaussian line were used to quantify the variation of the con- 
tinuum and Fe K emission over the remnant. Region L in Fig. || 
corresponds to the brightest hard continuum and a rather low Fe 
equivalent width and region P to rather weak continuum with 
the highest Fe K equivalent width. The same model was also 
used to fit the composite spectrum. The results are summarised 
in Table [l|. The apparent temperature in these fits is a mea- 
sure of hardness and is higher than the value reported above 
because the power law continuum has been omitted for sim- 
plicity. It is clear that the temperature for the composite fit is 
similar to that of the small regions (within the statistical accu- 
racy), again indicating that the hard flux is distributed through- 
out the remnant. As shown in Fig. ^| and quantified in Table |l| 
the variation in equivalent width of the iron line is very large. 
We attempted to fit the composite spectrum from MOS 1 + 2 + 
PN 4.07—1 2.0 ke V with a simple NEI model (see for example 
Borkowski 1994 ) plus a power law. The power law was again 
fixed at index 1.8 and the temperature kT fixed at 2.6 keV. This 
gave a good fit to the continuum with the Fe line blanked out. 

A reasonable fit was obtained for an ionisation timescale 
parameter of 5.0 x 10 11 s cm -3 or greater. For an age of 340 
years, (10 10 seconds) this implies a density of > 50 cm -3 
which is rather large. This fit gave about the correct equiva- 
lent width for the Fe line (assuming solar abundances), but did 
not fit the line profile very well, presumably because the line is 
Doppler shifted and/or broadened. It also gave a reasonable fit 
to the Ni line feature at around 7.7 keV. Lower more reason- 
able values of the time scale 5.0 x 10 10 s cm -3 corresponding 
to a density of 5 cm -3 predict a much larger equivalent width 
(a factor of 4 too large), a line energy which is too low and no 
Ni line at 7.7 keV. 

Using the same model on region P in the SE where the 
equivalent width of Fe is very high gave a poor fit. No com- 
bination of temperature or time scale could produce a strong 
enough line at the correct energy. The E knot is clearly anoma- 
lous with abnormally high Fe abundance and/or extreme non- 
thermal equilibrium ionisation. 
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3.4. Discussion 

The spectral form of the non-thermal high-energy "tail" is not 
inconsistent with a sim ple model of synchrotron emission from 
SNRs Reynolds 1998 . However in such a model the electrons 
need to be accelerated to energies of tens of TeV at the primary 
shock and the associated synchrotron emission would be ex- 
pected to be concentrated in the compressed magnetic field just 
inside the shock front. The hard X-ray continuum maps from 
XMM-Newton indicates that the 8.0—15.0 keV flux, predom- 
inantly due to the previously reported high-energy tail, does 
not originate from a few localized regions such as X-ray bright 
knots and filaments, nor does it originate from a limb bright- 
ened (fractionary) shell close to the shock front generated by 
the primary blast wave. In fact a low brightness outer ring 
structure, presumably associated with the primary shock, can 
be discerned in the hard X-ray image but it contains only a 
few percent of the total hard X-ray flux. Therefore the hard 
X-ray image is morphologically inconsistent with the simple 
synchrotron model developed by Reynolds ( 1998[ ). An alterna- 
tive explanation for the observation of hard X-ray tails in the 
spectra of supernova remnants is the presence of non-thermal 
Bremsstrahlung generat ed by a popul ation of suprathermal 
electrons (Asvarov et al. 1989| ; Laming 2000| ). One might ex- 
pect in this case some degree of correlation with the strong line 
emitting regions, however abundance variations make this po- 
tential correlation ambiguous. A search of Cas A for regions 
with spectra devoid of line emission, as a possible tracer of 
synchrotron emission, was unsuccesful although Chandra ob- 
servations indicated the presence of such a region just beyond 
the westernmost tip of the N rim of the remnant (Hughes et 
al. 2000 ). We specifically checked this region and still confirm 
the presence of, relatively weak, line emission, with for exam- 
ple equivalent widths of 200—300 eV for the He-like triplet of 
Si, S and Fe. We cannot exclude some contamination by line 
emission from neighbouring regions due to the wings of the 
XMM psf, on the other hand the Chandra data for this region 
also show line residuals above the continuum fit presented by 
Hughes et al. In summary we have not spotted any positive evi- 
dence, either morphologically or spectrally, for the presence of 
synchrotron emission in Cas A and, by implication, the associ- 
ated TeV electrons. Consequently, whether the hard X-ray tail 
is thermal or non-thermal remains an open question. 




5 10 15 

Wavelength (A) 

Fig. 4. First order (m = — 1) RGS spectra of Cas A for the three 
regions. Line identifications of the principle lines are provided. 





Fig. 5. RGS images of O VIII Ly-/3 (top left) and Ly-a (top 
right), the dispersion direction runs from top-left to lower-right 
at 45 degrees with the vertical. Emission is seen in three blobs 
in the E and N rim. The ratio map of the Ly-/? to Ly-o; im- 
ages is shown at the bottom left panel. The enhanced spot on 
the SE rim of the image indicates a Ly-a emission deficit in 
this region. For comparison we show a soft band MOS image 
(0.5 — 1 .2 keV) in the lower right panel. 



4. High-resolution spectroscopic data 

4. 1. Analysis & results 

Fig. [5| displays spectra from three extraction regions over the 
remnant situated in the N, NE and SE of the remnant. The 
data from RGS1 and RGS 2 have been combined in order to in- 
crease the statistical weight of the observation. Several lines of 
highly ionized species of Si, Mg, Ne, Fe L and O are detected in 
the spectrum. Due to interstellar absorption no features can be 
measured long-ward of ~20 A. The analysis and interpretation 
of these spectra will be the subject of a forthcoming paper. 

As a first result for this paper we extracted images of the 
O VIII Ly-a and Ly-/3 lines to probe small scale variations in 
absorption effects over this part of the remnant and to investi- 



gate the potential presence of resonance scattering in the limb 
brightened shells viewed edge-on. The temperature range rele- 
vant for Cas A has no influence on the Ly-a/Ly-/3 ratio. Since 
RGS is a slitless spectrometer it is possible to extract dispersed 
monochromatic images of Cas A. These images were con- 
verted from wavelength to spatial coordinates using the equa- 
tion AA = dsm(a)A(j), where Acj) is the offset along the dis- 
persion direction, AA the wavelength shift, a the angle of in- 
cidence on the gratings and d the line distance of the gratings. 
However, any Doppler broadening is also convolved along the 
dispersion direction and, depending on its magnitude, Doppler 
broadening could distort the RGS images. The O VIII Ly-a and 
Ly-/3 maps are presented in Fig. ||. The oxygen emission is seen 
to originate from three blobs on the E and N rim of the remnant. 
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Fig. 6. Upper panel: cross dispersion profiles of O VIII Ly-(3 
(thick line) and O VIII Ly-a (thin line). Lower panel: ratio of 
Ly-a to Ly-(3. Numbers indicate estimated AT H values from 
Keohane (1998), in units of 10 22 cm" 2 . 



The emission is faintest in the SE blob where there also seems 
to be a Ly-a deficit compared to Ly-/3. To investigate this we 
produced cross-dispersion image profiles of the Ly-a and Ly-/3 
lines by integrating over the dispersion direction. The cross- 
dispersion profiles together with the Ly-a to Ly-/? ratio have 
been plotted in Fig. || (note that this is the inverse of the ratio 
plotted in the brightness distribution in Fig. || because of the 
low Ly-a/Ly-/3 ratio). The ratio plot seems to steepen gradu- 
ally from a value of ~0.5 to ~0. 15—0.20, at a cross-dispersion 
above +1'. This decrease seems to indicate a Ly-a emission 
deficit in this region. The same effect is seen in the ratio image 
presented in Fig. ^. 

4.2. Discussion 

Three possible causes can be identified for explaining the 
O VIII Ly-a deficit in the SE part of Cas A, i.e. line blending, 
TVh variation and resonance scattering. 

(1) A line blending effect due to the reduced spectral reso- 
lution arising from the finite source extent causes the Ly-/3 line 
to be "polluted" by neighbouring lines from the Fe-L complex. 
Given the limited angular extent of the bright emission region 



in the dispersion direction, i.e. <1', and the fact that Cas A 
is an oxygen-rich remnant, this effect can be neglected, i.e. the 
oxygen emission lines strongly dominate the Fe-L lines and the 
spectral resolving power is still adequate (R~130). 

(2) The potential influence of variations in the column den- 
sity TVh in this region can be assessed by using the column 
densities given by Keohane (1998). These data take into ac- 
count both the H I and molecular absorption and have an av- 
erage resolution of 30". In the E part of Cas A they centre 
at 1.10 x 10 22 cm -2 with a full- width spread ranging from 
1.00 - 1.20 x 10 22 cm -2 (only the bright spot Cas A W shows 
a larger column of about 1.4 x 10 22 cm -2 ). Applying this full 
range going from N (least absorbed) to the SE bright rim (most 
absorbed) yields a decrease of the O vill Ly-a/Ly-(3 ratio of 
1 .4. The observed decrease of the Ly-a/Ly-(3 ratio, displayed in 
Fig. |6[ appears to be substantially larger, i.e. a factor >3 from 
about 0.45 to 0.15. Smaller-scale variations could be caused by 
small scale structure in the Nn value on a ~10" angular scale. 
To fully explain the decrease observed a column density varia- 
tion of 6 x 10 21 cm -2 (i.e. from 1.0 - 1.6 x 10 22 cm" 2 ) would 
have to be present, which seems rather unlikely. 

(3) Alternatively a decrease in the Ly-a/Ly- (3 ratio could 
also be introduced by the presence of resonance scattering of 
the Ovill Ly-a photons in the X-ray bright rims if viewed 
edge-on (Kaastra & Mewe 1995). A factor of 2— 3 reduction 
in Ly-a intensity can be explained by an optical depth in the 
Ly-a line of 2 or more (following the discussion in Kaastra & 
Mewe). For typical conditions in the rims of Cas A, an opti- 
cal depth larger than 2 is expected to be present if the micro 
turbulent velocity is smaller than ~700 km/s. 
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